Preprint typeset in JINST style - HYPER VERSION 



Study of the response and photon-counting 
resolution of silicon photomultipliers using a 
generic simulation framework 



Patrick Eckert , Rainer Stamen" and Hans-Christian Schultz-Coulon 

"Kirchhoff-Institut fiir Physik, Universitdt Heidelberg, 
Im Neuenheimer Feld 227, 69120 Heidelberg, Germany 
E-mail: Ipatrick . eckertBkip . uni-heidelberq . de| 



ABSTRACT: We present a generic framework for the simulation of Silicon Photomultipliers 
(SiPMs) which enables detailed modelling of the SiPM response using basic SiPM parameters 
and geometry as an input. Depending on the specified SiPM properties which can be determined 
from basic characterisation measurements, the simulation generates the signal charge and pulse 
shape for arbitrary incident light pulse distributions. The simulation has been validated in the 
whole dynamic range for a Hamamatsu S10362-11-100C MPPC and was used to study the effect 
of different noise sources like optical cross-talk and after-pulsing on the response curve and the 
photon-counting resolution. 
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1. Introduction 

Silicon Photomultipliers (SiPMs) are a modern type of solid state photo-detectors which are suit- 
able for a large range of applications such as calorimetry in high-energy physics, astrophysics or 
medical imaging [fll 0. ED- The key features of SiPMs are a very compact size (typ. 1 mm 2 ), insen- 
sitivity to magnetic fields, large gain, high photon detection efficiency and excellent single photon 
resolution. 

A detailed understanding of the response to photons and the dependence on the operation param- 
eters is essential for the utilisation of a photo-sensor. For SiPMs, the response depends on a mul- 
titude of effects like the photon detection efficiency, dark-rate, cross-talk, after-pulsing and pixel 
recovery time, as well as the spatial and time distribution of the incident light pulse. Due to the 
complex composition of the signal, an analytical description of the SiPM response is practically 
impossible. For this reason the Monte Carlo simulation framework GosSiP 1 was developed. It 
provides a detailed SiPM model which allows to predict the response under arbitrary operating 
conditions and offers the opportunity to study the interplay of the different effects influencing the 
SiPM signal; in particular, it enables to identify the signal component which limits the performance 
for a specific application. 

Earlier simulation studies presented in 0, ED already show a good description of the SiPM re- 
sponse for low light intensities. The simulation framework presented in this paper complements 
these studies by offering a more general simulation tool which is validated in the full dynamic range 
including the saturation regime at high light intensities. The new framework is used to study the 
effect of different noise contributions on the SiPM response and the photon-counting resolution. 

1 Generic framework for the simulation of Silicon Photomultipliers 
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2. SiPM Response Characteristics 

Silicon Photomultipliers consist of an array of typically several hundred up to several thousand 
micro-cells (pixels) per 1 mm 2 which are connected to a common output. Each pixel consists of 
an avalanche photo-diode (APD) operated at a bias voltage V^ius larger than the diode breakdown 
voltage Vtreak El- m this so-called Geiger mode operation, the electrical field at the p-n junction 
is large enough for free charge earners to initiate a self-sustaining avalanche process {avalanche 
breakdown). The avalanche is stopped with a large built-in quench resistor (passive quenching) 
or with an active electronic circuit controlling the bias voltage (active quenching). The presented 
SiPM simulation is designed for sensors with passive quenching, since this quenching scheme 
is implemented in most commercially available SiPMs. The following discussions are therefore 
restrained to SiPMs with passive quenching. 

A consequence of the Geiger mode operation is that the charge produced in the avalanche process 
is independent of the number of incident photons hitting a pixel. The photon-counting capability 
of SiPMs is hence achieved by segmenting the detector area into an array of individual pixels. The 
number of detected photons can then be determined from the number of pixels where an avalanche 
breakdown occurred. Besides the photon detection, there are however several other effects which 
contribute to the signal and thus have to be taken into account when analysing the SiPM response: 

• Gain: The gain of the SiPM determines the charge which is produced by a single avalanche. 
In good approximation, the gain depends linearly on the pixel capacitance C p i xe i and the 
applied over-voltage V over which is defined as the difference of the bias voltage and the 
breakdown voltage: 

Cpixd Vmer = _ (2 1} 

Here, q e is the elementary charge and C p i xe i the pixel capacitance. 

• Photon Detection Efficiency: An essential property of every photo-detector is the photon 
detection efficiency (PDE) which describes the probability to detect a single photon. For 
SiPMs, this is determined by the quantum efficiency, the geometrical fill factor, surface re- 
flections and the probability to initiate a Geiger discharge. The PDE significantly contributes 
to the photon-counting resolution due to the fluctuations in the number of detected photons. 

• Dynamic Range & Recovery Time: The finite number of pixels intrinsically limits the num- 
ber of photons which can be detected simultaneously. This leads to a saturation behaviour 
of the SiPM response and fundamentally limits the dynamic range of the SiPM as well as 
the photon-counting resolution for high light intensities. For an ideal sensor and an infinitely 
short light pulse, the SiPM response can be described as 

N av =N pixels -(l-e "W" ), (2.2) 

where N av is the number of avalanche breakdowns, N p i xe i s is the total number of pixels, Ny 
is the number if incident photons and Epde is the photon detection efficiency. Considering 
photons which arrive within a certain time interval, the pixel recovery time can significantly 
influence the SiPM response. During the avalanche and quenching process the bias voltage at 
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a pixel drops down to the breakdown voltage and recovers to the nominal operating voltage 
with a characteristic recovery time T recovery : 

Vover(t) = V over (0) • (1 - g-'/Wo-). (2.3) 

This recovery effect reduces the amount of charge which is produced in a subsequent avalanche 
breakdown. In first approximation 2 , the recovery time is determined by the pixel capacitance 
and the quench resistor: 

^recovery — Rquench ' C pixel ■ (2-4) 

• Cross-talk: It is well known that charge carriers crossing a p-n junction with reverse bias 
voltage beyond breakdown can create photons in the visible range fl7[|. The probability to 
emit a photon with sufficient energy to create an electron-hole pair is about 3 x 10~ 5 per 
charge carrier crossing the junction. For a SiPM with a typical gain of 10 6 about 30 photons 
are produced per avalanche. These photons can propagate into nearby pixels and initiate a 
secondary avalanche breakdown faking a "real" photon hit. This effect is referred to as opti- 
cal cross-talk and degrades the photon-counting resolution of a SiPM due to the fluctuations 
in the number of cross-talk events. In addition, the increased pixel occupancy caused by 
optical cross-talk reduces the dynamic range of the sensor. 

• After-pulses: During an avalanche, charge earners can be trapped in the silicon lattice and 
be released after a characteristic time. This can cause a delayed second avalanche breakdown 
faking a photon signal and thus degrading the photon-counting resolution. In principle there 
can be several kinds of trapping centres with different characteristic trapping times. The data 
shown in section ^TT] and ED indicate the existence of two kinds of trapping centres, one with 
a slow and one with a fast time constant. Since cross-talk and after-pulse events are both 
triggered by an avalanche breakdown, these effects are often labelled as correlated noise. 

• Dark-rate: Avalanche breakdowns cannot only be triggered by photon absorption but also 
by thermal or field mediated excitations of electrons in the silicon lattice. These events will 
in the following be denoted as thermal pulses. The thermal pulses and secondary avalanches 
caused by cross-talk and after-pulses make up the dark-rate of a sensor. 

• Electronic Noise & Excess Noise: The single photon resolution is determined by the pedestal 
noise o ped and gain fluctuations 

"gain ■ 

o( i ) 2 = Oped + i -Og ain , (2.5) 

where i is the number of avalanche breakdowns. The pedestal noise is defined by the signal 
fluctuation when no avalanche breakdown occurs. These fluctuations arise from leakage 
current and noise from the readout electronics. The gain fluctuations originate from pixel-to- 
pixel variations of the gain due to variations in the quenching resistance and fluctuations in 
the avalanche process. The value o ga j n /G is often referred to as Excess Noise Factor. In the 
following, a p e d and o ga i„ will be denoted as Electronic Noise and Excess Noise, respectively. 

2 Valid for small input impedance of the readout electronics 
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3. SiPM Simulation Framework 

The prediction of the SiPM response for arbitrary operation conditions is non-trivial, due to the 
complex composition of the SiPM response - especially in the case of large correlated noise, high 
light intensities and arbitrary light pulse shapes. For this reason, the Monte Carlo simulation frame- 
work GosSiP 3 , which provides a detailed model of the SiPM response, has been developed. The 
simulation can be adjusted to any type of sensor by providing the basic SiPM parameters and pixel 
arrangements. In the following sections, this simulation framework and its comparison to data is 
presented. 

The simulation framework is divided into three main parts and a graphical user interface which 
allows to adjust the SiPM and light source parameters, select a measurement and display the simula- 
tion result. Due to the modular structure the integration of GosSiP into other software frameworks, 
like e. g. Geant4 4 , is straightforward. 

• Light source: A basic simulation of a light source is implemented which generates a cus- 
tomisable light pulse illuminating the SiPM. The adjustable parameters are (i) the number 
of photons with the option of a Poisson smearing, (ii) the duration and time structure of the 
pulse (flat, Gaussian or exponential decay) and (iii) the position and shape (rectangular or 
elliptic) of the light spot. The position and time stamp of the individual photons within the 
pulse are randomly distributed according to the specified time and spatial structure. 

• SiPM response: The centrepiece of the GosSiP framework is the simulation of the signal 
charge and waveform for the specified SiPM parameters and light pulse input. All SiPM 
parameters can be customised in order to model different types of SiPMs. 

• Data acquisition: A simple simulation of the basic characterisation measurements using 
standard electronics is provided, including measurements of the charge and time spectra, 
dark-rate threshold scans and response curve determination. 

In the first simulation step, a list of incoming photons is generated by the internal light source 
simulation. Alternatively, this list can also be generated by external simulation tools. The photon 
information is then passed to the SiPM simulation which models the response of the sensor to 
the incident light pulse. Each photon hitting the sensor can trigger an avalanche with a certain 
probability given by the PDE. To simulate the thermal noise, random pulses are generated according 
to the specified thermal pulse time constant (see equation pT71 ). The information about the initial 
time of each avalanche as well as the position of the fired pixel is stored in a list which is then 
processed chronologically according to the following steps. 

The charge generated by the avalanche is determined using the specified gain and a Gaussian 
smearing to account for the excess noise. Then cross-talk and after-pulse events are generated 
according to the specified cross-talk and after-pulse probabilities (see equation |7B] & |4.10D and 
inserted chronologically into the avalanche list. This naturally enables higher-order noise cascades 
from an initial avalanche. A simple cross-talk model is used in which cross-talk can only occur in 

3 www.kip.uni-heidelberg.de/hep-detektoren/gossip 
4 http://geant4.cern.ch 
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the four neighbouring pixels. This model yields a good description of the data as shown in earlier 
publications 0, @H and validated below. The after-pulse model is based on the probability density 
function given by equation which depends on a characteristic trapping time constant and a trap- 
ping probability. Two independent after-pulse components are implemented accounting for a slow 
and a fast trapping time constant. 

The chronological processing of the avalanche list is essential in order to model the pixel recov- 
ery effects. The recovery process influences all effects which contribute to the signal (i. e. gain, 
excess noise, PDE, thermal rate, cross-talk, after-pulses) due to the changing over-voltage (see 
equation |2~3"1 ). The recovery of the gain is modelled using equation ^3] and [2~T1 : 

G{t) = G(0) • (1 -g-'A-'—y). (3.1) 

The recovery of the PDE, excess noise, cross-talk, after-pulse and thermal pulse probabilities is 
modelled in the same way. This is only a rough approximation since the latter effects in general 
do not depend linearly on the over-voltage, as it is the case for the gain. However, the impact of 
this approximation is small due to the reduced signal amplitude of a firing pixel which is not fully 
recovered. A more sophisticated recovery model can be implemented easily. This would how- 
ever necessitate a more complex simulation input, since the over-voltage dependence of all effects 
would have to be provided. 

In the final step of the simulation chain, the simulated output signal can be processed by the sim- 
ulation of the data acquisition which models the most common characterisation measurements. 
Alternatively, the simulated output can be passed to an external software for a custom digitisation 
of the signal. 

4. Simulation Input & Validation 
4.1 Input Parameter Determination 

The simulation input is given by the following basic SiPM parameters: PDE, gain, excess noise, 
electronic noise, thermal noise time constant, cross-talk probability, after-pulse probabilities and 
after-pulse time constants, recovery time, number of pixels, the pixel arrangement and the single 
pixel pulse shape. These parameters can be determined from basic characterisation measurements 
which will be briefly described in this section. A detailed description of the measurement setup 
and method can be found in [JJ. The following study was earned out using a Hamamatsu MPPC 5 
S 10362- 11-100C device with 100 pixels, operated at an over- voltage of 0.5 V and 1.0 V. This covers 
the range of low and high noise operation of the sensor. Table [I] summarises the obtained values 
for the input parameters. 

Quench Resistance & Recovery Time The quench resistance R quenc h is determined from a linear 
fit of the I-V curve for forward bias operation (see figure [I]): 

R = Rquench/ N pUe l s + R C = AV/AZ, (4. 1) 

5 Multi Pixel Photon Counter 



-5- 



Table 1. Simulation input parameters determined by the characterisation measurements. The PDE is mea- 
sured at X s» 658 nm (red laser diode), whereas the maximum sensitivity of the used sensor is at A « 460 nm 
(see Ml 



Parameter 


V over = 0.5V 


V over = 1.0V 


PDE [%} 


4.8±0.3 


11.6±0.6 


Gain [q e ■ 10 5 ] 


9.2±0.2 


21.5±0.2 


Ogain/G&m [%] 


13.2±0.3 


5.8±0.3 


a ped /Gam[%] 


6.4±0.3 


9.1 ±0.3 


T DR [ns] 


6632 ±66 


2821 ±28 


T AP t , ol ,, H 


211 ±26 


181 ±6 


*AP fmt H 


61±7 


61±2 


[%] 


2.3±0.1 


13.0±0.6 


PAP f « [%} 


3.6±0.2 


16.3±0.4 


P C T,n>l [%] 


2.4±0.1 


14.3±0.6 


^recovery [is] 




38.4±1.2 



where Roc is a serial resistance in the operation circuit. For the tested sensor a value of R que nch = 
96 ± 3k£2 is obtained. This yields an indirect measurement of the pixel recovery time using equa- 
tion ^3] and a pixel capacitance of C p i xe i = 400 ± 1 fF which is determined from the gain - voltage 
dependence (see section |Q]) The resulting recovery time is X recovery = 38.4 ± 1.2ns which is com- 
parable to the value of ~ 33 ns presented in a study UTOl based on a waveform analysis. 

x10" 3 

0.4 
0.35 
0.3 

< 0.25 
S 0.2 
= 0.15 
0.1 
0.05 


1 2 3 4 5 

Forward Bias Voltage [V] 

Figure 1. IV-Curve: The quench resistance determined from a linear fit to the the data above 1 V forward 
bias voltage as in this region the resistance of the pn-junction is negligible. 

Charge Spectrum The charge spectrum is one of the most important tools for the characteri- 
sation of SiPMs, since it provides information about a multitude of parameters. It is acquired by 
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Figure 2. Single photo-electron spectrum for 0.5 V over- voltage. Each peak corresponds to a certain number 
of photo-electrons (pe). 



illuminating the sensor with a short light pulse (^?(ns)) and integrating the signal charge within a 
certain time window using a QDC 6 . A typical charge spectrum for low light levels (single pixel 
spectrum), as displayed in figure 0, shows several well separated peaks which correspond to a cer- 
tain number of avalanche breakdowns. All events with no pixel signals are contained in the first 
peak (pedestal), the second peak contains all events with one pixel firing, and so on. The SiPM gain 
is given by the charge which is produced by a single avalanche and hence is equivalent to the charge 
difference between two peaks. The voltage dependence of the gain described by equation ^71] can 
be used to determine the pixel capacitance and breakdown voltage. 

The number of photons in the incident light pulse, and therefore also the number of detected pho- 
tons, is expected to be Poisson distributed. However, the distribution of the number of detected 
photons is distorted due to cross-talk and after-pulses which fake photo-electron signals. The mean 
value of the charge spectrum therefore corresponds to the mean number of avalanche breakdowns 
N av including signals from cross-talk, after-pulses and thermal pulses 

Nav = (Q-Q P ed)/G, (4.2) 

with Q denoting the mean value of the charge spectrum and Q pe d denoting the pedestal value. 
Exploiting Poisson statistics, the mean number of detected photons N pe without any bias from 
noise can be determined from the number of pedestal events which are not affected by cross-talk 
and after-pulses: 

P(0,N pe ) = e- N ?° (4.3) 
N pe = -ln(P(0,N pe )) = -ln( CDR M Nped ). (4.4) 

6 Charge-to-Digital converter; LeCroy Model 2249, 50 £1 input impedance. 
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Here, N tot is the total number of recorded events and N pe d is the number of pedestal events, c^r is 
a correction factor rescaling the value of N pe d in order to account for the dark-rate which reduce 
the number of pedestal events. This correction factor can be determined from a charge spectrum of 
the dark-rate [0]: 



where cdr — 1 corresponds to the probability for a dark-rate event to occur during the charge inte- 
gration. 

If the charge spectrum is measured using a referenced light source, this statistical analysis allows 
to determine the PDE without the effects of cross-talk, after-pulses and dark-rate. A detailed de- 
scription of this measurement is given in [0). 

If no pixel fires only the charge generated by electronic noise is integrated; this is reflected in the 
width of the pedestal peak o pe d- The excess noise is given by the additional broadening of the 1 pe 
peak compared to the pedestal: 



where 0\ pe is the width of the 1 pe peak. 

Dark-rate Time Spectrum The thermal pulse and after-pulse characteristics can be determined 
from a statistical analysis of the time interval between two consecutive dark-rate pulses. The spec- 
trum of this time interval is acquired using a discriminator at 0.5 pe level and measuring the time 
stamps of the discriminator pulses with a TDC 7 module. For our setup, this method allows to stably 
acquire time intervals down to t « 100ns, due to the relatively wide pulse shape of the signal. The 
time spectrum represents the probability density for a dark-rate pulse to occur after a time t with 
respect to the previous pulse. The probability density for the thermal pulse component is given by: 

1 __L_ 

PTp{t) = e 'tp, (4.7) 

Z T p 

where Xjp is the characteristic thermal pulse time constant. For after-pulses with a time constant 
of Tap, this probability is multiplied with the after-pulse probability Pap'. 

P ap t_ 

p A p{t) = — -e *ap. (4.8) 
Tap 

The measured time spectrum can be approximated with a combination of the thermal pulse, slow 
and fast after-pulse component: 



c DR = N™/Ni 



(4.5) 




(4.6) 



N(t) =N- 



(4.9) 




+ PAP f (t)-[l 



f 

Jo 



p T p{t')dt'}-[l 



f 

Jo 



PapM') dt'] 



) 



7 Time-to-Digital Converter; CAEN V1290A 
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Figure 3. Spectrum of the time intervals between two consecutive dark-rate pulses. 

where pAp, f (t) denotes the probability density for the slow and fast after-pulse component, respec- 
tively. The simulation input parameters XTp,lAP s ,lAP f ,PAP s ^^PAP f can be extracted by fitting the 
measured time spectrum with formula However, formula is a first order approximation and 
does not take into account cross-talk and pixel recovery which influence the number of after-pulse 
events. Therefore, the values for the time constants and after-pulse probability determined by the fit 
exhibit a systematic uncertainty. The minimum time interval of t « 100 ns which can be achieved 
with our setup also limits the sensitivity to fast after-pulses which may lead to systematic errors 
in the fit. In order to account for these effects, the SiPM simulation output is fitted to the data by 
iterative variation of the input parameters. The resulting values for the input parameters show a 
deviation of up to 20% compared to the values obtained from the analytical fit using equation P~9l 
With a more sophisticated setup for the determination of the after-pulse properties (e.g. [2J, 0]) 
it should be possible to determine the simulation input parameters via an analytical fit similar to 
equation 

Dark-rate Threshold Scan The cross-talk probability can be determined by counting the number 
of dark-rate events with a certain number of fired pixels. This is done using a discriminator module 
which selects dark-rate pulses with a certain amplitude. The discriminator signals are then counted 
with a scaler module 8 . Figure § shows the dark-rate as a function of the applied discrimination 
threshold. The total dark-rate at 0.5 pe threshold is composed of the thermal pulse rate and a 
contribution from after-pulses. The events where two or more pixels fire at the same time, which 
correspond to a discrimination threshold of 1.5 pe, are caused by dark-rate events triggering optical 
cross-talk. The probability for two thermal pulses to occur at the same time is negligible. The 
probability for an avalanche to generate one or more cross-talk events PcT,n> l is given by the ratio 
of the dark-rates at 1.5 pe and 0.5 pe threshold: Pcr,n>i = DRi.5 pe /DRo.5 pe ; this is the cross-talk 

8 LeCroy Model 2250B 
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Figure 4. Dark-rate measured at different discrimination thresholds. 



input parameter for the simulation. A simplified cross-talk model is implemented in the simulation 
taking into account only the four directly neighbouring pixels which generate a cross-talk event 
with a respective probability of Per- This value is related to the input value via the following 
equation: 

(l-/ 5 cr) 4 = l-Pcr ! „>i. (4.10) 

Simulation Cross-checks Comparing the simulated response to the characterisation measure- 
ments is important to validate the physics models used in the simulation. The shape of the simu- 
lated time spectrum (see figure 0) agrees well with the measurement confirming the implemented 
after-pulse and thermal pulse model. The simulation also reproduces the total dark-rate at 0.5 pe 
threshold shown in figure |] which independently validates the after-pulse and thermal pulse model. 
The rates at 2.5 pe and 3.5 pe threshold are slightly over-estimated by the simulation; this could 
be attributed to the simplified cross-talk model or the simplified discriminator model 9 used in the 
simulation. However, the deviation is small and the probability PcT,n>2 for two or more cross-talk 
events which corresponds to the rate at 2.5 pe threshold is smaller than 2%. Therefore, the effect 
can be safely neglected. 

The simulated charge spectrum (see figure 0) also agrees well with the data confirming the im- 
plemented gain, excess noise, and photon detection model. Since the amplitudes of the individual 
peaks in the spectrum are also influenced by cross-talk and after-pulses, this again shows the valid- 
ity of the implementation of these two effects. 

These cross-checks verify, that the simulation offers a good description of the SiPM dark-rate and 
the response to low light intensities. Similar results with a comparable simulation approach have 
been published in 0]. However, for large light intensities, the high pixel occupancy and pixel 
recovery have a significant influence on the response and effects like cross-talk, after-pulses and 

9 The discriminator is modelled using a simple parametrization. 
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Figure 5. Schematic view of the setup for the measurement of the response curve. 

photo detection. Therefore, the validation of the simulation was extended to cover the whole dy- 
namic range. This will be discussed in the following chapter. 

4.2 Response Curve & Resolution 

For low light intensities, the individual pixel signals can be treated as independent which allows a 
rather simple description of the SiPM response. This is not the case for high light intensities due 
to the saturation of the signal. In order to validate the simulation over the whole dynamic range 
including the low intensity region as well as the non-linear saturation region, the response curve 
and signal resolution were measured and compared to the simulation results. The response curve 
is a sensitive quantity, since it is influenced by all SiPM parameters which enter the simulation. 

The basic principle of the measurement is to determine the SiPM response as a function of 
the incident photon flux (see figure |5|). The SiPM is illuminated using a laser diode emitting at a 
wavelength of A « 658 nm, which is driven by a pulse generator with 4 ns pulse width. A beam 
splitter evenly distributes the light to the SiPM and a NIST 10 calibrated photo-diode which serves 
as a reference sensor with ideal linear response. The charge of the SiPM signal is measured with a 
QDC within an integration gate of 300 ns. 

The light spot properties are crucial for the measurement of the response curve. In order 
to achieve a homogeneous illumination of the SiPM, a diffuser is used to obtain a uniform light 

10 National Institute of Standards and Technology 
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Figure 6. Leading edge of the signal (about 20 firing pixels) for data (red), reconstructed waveform assuming 
a flat 1.0 ± 0.5 ns light pulse (black line and green uncertainty band), and reconstructed waveform assuming 
a fiat 4.0ns light pulse (black dashed line). 



spot. In addition, the MPPC S10362-1 1-100C sensor was chosen for this measurement due to 
its long pixel recovery time of 38.4ns which minimizes the influence of the time structure of the 
light pulse. However, it was observed, that the SiPM response in the non-linear range still depends 
weakly on the width of the voltage pulse driving the laser diode. Since the time response of the laser 
diode to the voltage pulse is not known, the duration of the light pulse was estimated by analysing 
the leading edge of the SiPM signal at medium light intensities (« 20 firing pixels). The signal 
waveform was reconstructed using the simulation adding up averaged 1 pe waveforms which have 
been measured with an oscilloscope. The best agreement of measured and reconstructed signal 
was achieved with a light pulse time interval of At « 1.0 ± 0.5 ns and a flat time distribution (see 
figure ^). Assuming a pulse width of 4 ns, which corresponds to the pulse width of the pulse 
generator, a significant deviation between measured and reconstructed signal is observed. 

The number of photons, N r , hitting the sensor is determined with the NIST reference sensor. 
The sensor is calibrated using SiPM charge spectra at low intensity, from which N r is obtained 
via the relation N r = N pe /PDE. N pe , the average number of photons detected with the SiPM, is 
extracted via equation and the PDE is taken from Table [I]. The uncertainty of this calibration 
procedure leads to the systematic uncertainty on the obtained value of Ny. 

The SiPM response, i. e. the number of fired pixels Nov, is determined from the measured signal 
charge using equation Due to the limited dynamic range of the QDC the response at high 
light intensities cannot be measured using pre-amplification. The necessary gain factor (in units of 
QDC-channels), G, for non-amplified signals thus has to be determined. This is done by comparing 
the measured charge Q — Q pec i at low light levels with and without amplification. The ratio of the 
measured charges yields the amplification factor which is used to correct the gain measured using 
amplified single photon spectra. The uncertainty on this amplification factor leads to the dominant 
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systematic error on N av . 

Figures [7j PL R [TD| show the simulated and measured response curve and the resolution for 
0.5 V and 1.0 V over-voltage. The contribution from the dominant noise sources determined by 
the simulation are shown beneath. The error band of the simulated response originates from the 
measurement uncertainties of the input parameters. 

For 0.5 V over- voltage, the response curve is fitted using a formula similar to equation [ 



C-Ny 

N av {N y ) =N pixels • (4.11) 

where e describes the mean number of avalanches, including PDE, cross-talk and after-pulse ef- 
fects, which are triggered by an incident photon. Ndr is the mean number of dark-rate events 
occurring within the time interval of the integration gate. For 1 .0 V over- voltage, this formula has 
to be modified due to the large contribution from cross-talk and after-pulses UTTH : 

N av (N 7 ) =N pixeh -{\-e~^)/{\-e C N-e~^)+N DR , (4.12) 

where Scn describes the contribution from correlated noise. 

For both low and high over-voltage, the simulation reproduces the measured response and 
resolution within the uncertainties and thus validates the simulation with an accuracy of about 5%. 
For the specific SiPM used (Hamamatsu MPPC S10362-11-100C) the simulation identifies after- 
pulsing as the dominant source of noise. The contribution from cross-talk to the SiPM response 
and resolution in the non-linear range is suppressed due to the high pixel occupancy. 

The relative photon-counting resolution On /Ny of a SiPM is determined by the inverted re- 
sponse curve Ny(N av ) and the resolution a^ m : 



On 7 _ dNy/dNw ■ <T Nm 



Ny Ny 



(4.13) 



The relative photon-counting resolution determined via this formula is shown in figure |TT] and [j^. 
Assuming no saturation effects, the relative resolution can be described by 

°Ny A B fA1AS 

= © , (4.14) 

Ny Ny ^JNy 

where A contains the contributions from the dark-rate, and B represents the stochastic term includ- 
ing fluctuations of the number of detected photons due to a limited PDE as well as contributions 
from cross-talk and after-pulses. Figures |TT] and [12] also show a fit of this function to the measured 
resolution using only the low intensity data points corresponding to a response below 20 pixels. 
The fit is then extrapolated to high intensities. It can be seen, that for the tested sensor, this descrip- 
tion is valid up to « 40 firing pixels, which corresponds to « 1000 photons for 0.5 V over-voltage 
and « 300 photons for 1.0 V over- voltage. For higher intensities, the signal saturation contributes 
significantly to the relative resolution. In case of 0.5 V over-voltage, the sensor performs well in the 
whole measurement range up to 4000 incident photons while the best relative resolution of 15% is 
achieved at « 2500 photons 70 fired pixels). For 1.0V over-voltage, a similar value of 20% at 
« 80 fired pixels corresponding to « 1000 incident photons is achieved. For > 100 firing pixels 
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the relative resolution is degraded due to the strong saturation of the signal. 
The simulation shows, that in a wide range the dominant contribution to the relative resolution is 
the limited PDE, whereas the effects of cross-talk and after-pulsing are small. Only in the saturation 
region cross-talk and after-pulsing contribute significantly to the relative photon-counting resolu- 
tion, due to the large fluctuations in the number of after-pulse and cross-talk events. In addition, 
the cross-talk also reduces the dynamic range due to the increased pixel occupancy, which further 
degrades the relative resolution in the non-linear range. 

5. Summary 

In this paper, a simulation framework of the response of Silicon Photomultipliers is presented. 
The tool offers a generic and detailed modelling of SiPMs and allows easy integration into custom 
applications. The simulation is validated in the full dynamic range for a 100 pixel Hamamatsu 
MPPC sensor operated at low and high bias voltage. The data are described by the simulation within 
approximately 5% uncertainty. The photon-counting resolution of the sensor is measured as a 
function of the light intensity and the contributions from cross-talk and after-pulses are determined 
using the simulation. 
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Figure 10. Top: Resolution of the response for high noise (1.0 V over- voltage). Bottom: Contribution from 
the individual noise sources. 
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Figure 11. Top: Photon-counting resolution for low noise (0.5 V over- voltage). The dotted line shows the 
resolution assuming no saturation effects. Bottom: Contribution from the individual noise sources. 
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Figure 12. Top: Photon-counting resolution for high noise (1.0 V over- voltage). The dotted line shows the 
resolution assuming no saturation effects. Bottom: Contribution from the individual noise sources. 
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